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Abstract. Comparison between the exact value of the spectral zeta function, 
Zh{1) = 5-^/^ [3 - 2cos(7r/5)] r2(l/5)/r(3/5), and the results of numeric and WKB 
calculations supports the conjecture by Bessis that all the eigenvalues of this PT- 
invariant hamiltonian are real. For one-dimensional Schrodinger operators with 
complex potentials having a monotonic imaginary part, the eigenfunctions (and the 
imaginary parts of their logarithmic derivatives) have no real zeros. PACS numbers: 

03.65.Db, 02.30.Tb, 02.60.Lj, 11.30.Er 



1. Introduction 

Some years ago Bessis (1995) had conjectured that all the eigenvalues of the cubic 
oscillator with purely imaginary coupling constant, 

H^^^+igx^; lni(g) = 0, (1) 

were real. Numerical and WKB investigations have found only real eigenvalues (Bessis, 
1995, Bender and Boettcher 1998ab, Bender et al 1999), although a rigorous proof is 
still lacking. 

The operator (0), defined on the intersection of the domains of d^/dx^ and 
is closed and has a compact, Hilbert-Schmidt resolvent (Caliceti et al 1980, Caliceti 
1999). Thus, its spectrum consists only of eigenvalues, which can acumulate only near 
infinity. Since its numerical range lies in the right complex half-plane [Re(z) > 0], the 
real parts of the eigenvalues are positive (Kato 1966). 

A Schrodinger operator 

H = -^ + Uix), (2) 

is called PT-invariant if the potential satisfies 

U{x) = U{-x). (3) 

Here, and in the following, we use the overline to denote complex conjugation. 
Although (|l|) is not Hermitean (symmetric), it is PT-invariant, but not separately 
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invariant under the spatial inversion P, [(t>{x) — > (j){—x)], or the time reversal T, 
[4>{x) (pix)], which map it on its adjoint H* . 

Adding harmonic term to the potential in we obtain 

U{x) = igx^ + iJLx^\ A* > 0. (4) 

The PT-invariant operator defined by (Q) and (^) has been intensely investigated. 
Caliceti et al (1980) have shown that each eigenvalue of the harmonic oscillator (to 
which it reduces for g = 0) stays real for sufficiently small real g or, equivalently, for 
sufficiently large (Delabaere and Pham 1999). 

The resolvent [z — H)^^ of the one-dimensional Schrdinger operator (^) is an 
integral operator for all z in the resolvent set. Its kernel, Gz{x,y) can be expressed 
through 2 solutions that satisfy the boundary condition at only one end of the interval 
(Naimark 1964, Caliceti et al 1980). Let f^{x;z) [ f-{x;z) ] be the (unique, up to a 
multiplicative constant) solution of the equation 

(z - H)f = f" + [z~ U{x)]f = 0, (5) 

which, together with its first derivative, is locally near +oo (respectively —oo). 
One can readily check that the kernel is 

r, f \ S{x-y)f+{x;z)f_{y,z)+0{y-x)f+{y;z)f^{x;z) 

""'^^''^^ w[.u,UKz) ■ 

Here 6{x) is the Heaviside function (characteristic function of the positive half-axis), 
and 

W[f+J-]iz) = f[{x; z)f^{x; z) - fLix; z)f+{x; z), (7) 

is the (x-independent) Wronskian of and /_. 

Generally speaking, a PT-invariant Schrodinger operator (|^), (|^) can have any 
kind of spectrum, as long as it is invariant with respect to complex conjugation. For 
example, let 

Ug,a — i5tanh(a;/a), with g, a> 0. (8) 

This PT-invariant potential of can be thought as an interpolation between the 
potential iFx of an imaginary electric field (g = aF, a — > -l-oo) and the imaginary step 
potential igsign(a;) (in the limit a ^ 0). Define Hg a by (^) with the potential (||) on 
the domain of the second derivative. The solutions of the eigenvalue equation for (||) 
can be obtained by analytic continuation from those given in (Landau and Lifschitz 
1977). 

The operator Hg a has no eigenvalues. Its spectrum, which does not depend on 
a, is the union of two half-lines parallel to the real axis: {±ig + 1; t > 0}. In the 
imaginary electric field limit, the lines move to infinity, so that the spectrum of the 
limiting operator reduces to the point at infinity (see e.g. Herbst 1979). 

Another instructive exactly soluble model is the PT-invariant "complexified 
harmonic oscillator" with the potential 

r &^x^, if a; > ; , . 

= \ e-'^x^ ifx<0; 

where |A| < 7r/2. For A = 0, (^ coincides the potential of the harmonic oscillator, 
while for A — ±7r/2 it is purely imaginary and antisymmetric, mimicking the cubic 
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potential in Let us also note that, in the limit A ^ the perturbation (^) on 
the harmonic oscillator is somewhat milder than for the e — > limit of the potential 
— (ia;)^+'^ « a;^[l + eln(|a;|) + i^sign(a;), considered by Bender and Boettcher(1998a). 

The operator H\, defined on Z)(d^/da:^) n D{x^), has compact (Hilbert-Schmidt) 
resolvent. But, for any A ^ only a finite number of the eigenvalues are real. The 
rest come in complex-conjugated pairs. 

In the next section, we will compute exactly one value of the spectral zeta 
function the sum of the inverse eigenvalues, and will compare the result 

with numerical and WKB calculations for the real eigenvalues of (|l|), obtaining a 
remarkable coincidence. Since the real parts of any presumed complex eigenvalues 
must be positive, this places stringent bounds on the values of the maximum value of 
the ratio \Re{Ea) / E^] for complex eigenvalues. In section 3 we will also show that 
the eigenfunctions (and their first derivatives) have no zeros on the real axis. This 
property holds for all potentials with monotonic imaginary part. 



2. Sum of inverse eigenvalues 



In our case we can set g — 1 without altering generality. The eigenvalue equation (^ 
reads 

/" +{z~ ix^) / = 0. (10) 

Generally speaking, the solutions of (|lO|) cannot be expressed in terms of the classical 
special functions. This is possible only for z ~ (Gradshtein and Ryzhik 1965, Caliceti 
et al 1980). Then, we can readily verify that 



and 
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X <0: 



(11) 



f^{x;0) ^ f+i-x;0), (12) 

are a pair of PT- conjugated solutions of the Equation ([l0|), and that f+{x; 0) [f-{x; 0)] 
and their first derivatives are locally near +00 [—00]. Here, I^, and K^, denote the 
Bessel function of the first kind and imaginary argument, respectively Macdonald's 
function. 

Although this is not immediately apparent from the representation (|l^), the 
solutions /+(x;0) and /_(a;;0) are integer functions of x. The first terms of the 
McLaurin series for /+(x; 0) are 
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^x + 0{x^) 



From ( p^ ) and (12) we get the z — value of the Wronskian (|^) 

w[f+j-m 
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The asymptotic behaviour of f+{x; 0) near ±00 is given by 



V5tt 
2x4. 



exp ( 



|eTa:5 - i^) 



) l + 0{x-i) 



exp(|e-x|a;|i + ^) l + 0{\x\-i) 



X 
X 



+00; 
—00. 



(13) 



(14) 



(15) 



Using the asymptotic behaviour (p^), Cahceti et al (1980) have shown that the 
kernel Go{x,y) belongs to the Hilbert-Schmidt class. 



|Go||hs 



oo r+oQ 



da;dy|Go(x,y)p < cxo. 



(16) 



Here we will note that the diagonal elements of the continuous in {x, y) Hilbert-Schmidt 
kernel Go{x,y) is absolutely summable, 



H-oo 

dx\Go{x, x)\ < oo, 

-oo 

since the continuous integrand is summable near both infinities, 
\f+{x;0)f-{x; 0)1 « O (la^p^/^^ , whenx ^ ±oo. 
The trace is equal to the sum of all the inverse eigenvalues 
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(17) 



(18) 
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which is finite. Analogous formulae can be written for the zeta functions of integer 
arguments s > 1. Taking into account the PT-conjugacy of /+ and we obtain 



. . 2sin^ 
Zh{1) - ^^Re 



AxU{x;0)f.{x-Q) 



4-00 



56/5r(f) 



dt t~^Ki{t) 



3 — 2 cos ■ 



cos(7r/10) = 
2.835094933. 



(20) 



In the last integral we rotated the contour of integration by — tt/IO, making the 
substitution x = (5t/2)'^ e""^/^", which leads to tabulated integrals (Gradshtein and 
Ryzhik 1965, 6.576.5). 

Bender and Boettcher (1998a) gave a WKB estimate for the large eigenvalues of 

® 



5, - 
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Comparison of (|T]) with the numerical results of Bender and Boettcher (1998a) shows 
that even for j = the error is about 5%. 

The sum of all inverse WKB eigenvalues ( pT| ) can be expressed in terms of the 
Riemann zeta function: 



'^'''=^11 -V—— - 



-r(|) 



3r(i) 



6/5 



2.885673793. 



(22) 



We can correct the WKB approximation ( p^ ) for Zh{^) by taking the numerical 
values of the first five eigenvalues from Table 1 of Bender and Boettcher (1998a), 
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who give their result to five digit precision. This yields 2.8351, to be compared with 
2.835094933 from (20). A remarkable coincidence! The two values coincide within the 
five-digit acuracy of the numerical eigenvalues. Since the real parts of any presumed 
complex eigenvalues must be positive, no compensations can occur in the sum. If some 
complex conjugate eigenvalues Ea are present, then the largest ratio 

Re(Sa) 



max 



(23) 



We consider this as strongly supportive for the conjectured (Bessis 1995) absence of 
complex eigenvalues. 



3. Nodeless character of the eigenfunctions 

As we can see from (^), the sum ( |l9| ) converges slowly. An improvement in the 
precision of the comparison between (20) and numerical result swill require precision 
computations of many eigenvalues. Bender and Boettcher (1998a) integrated the 
differential equation along the anti-Stokes lines in the lower complex half-plane, 
starting with large values of \x\, where the asymptotic behavior is attained, and 
matching the solutions at a; = 0. The following proposition on the properties of 
the eigenfunctions could prove useful for simplifying numerical work. 

Let (pE he an eigenfunction corresponding to the eigenvalue z ~ E in ([l0|). Then, 
proceeding as in the case of proving current conservation for the solutions of the 
stationary Schrodinger equation, we obtain 



j'e{x) = Im [cl>'E{x)<j>E{x)\ = [x' - lm{E)]\<j>E{x)\'. (24) 

The current density associated to an eigenfunction is a decreasing (increasing) function 
of a; if a;"^ — Im{E) < (>)0. Since the limit j£;(±oo) = 0, the current density is strictly 
negative on (— oo, +oo) and has a unique minimum at a; = y'hn{E). Thus, 

a. The eigenfunctions of (\1(\) have no zeros on the real axis; 



b. The imaginary part of their logarithmic derivatives, lm[(f)'^{x)/(f)E{x)], have 
constant sign. 

This proposition can be readily generalized to the eigenfunctions of Equation ^ with 
a complex potential U{x) having a monotonic imaginary part. 

Using the above proposition we can look for the eigenvalues of (^0|) by a shooting 
method, integrating the Ricatti equation satisfied by the logarithmic derivative 
s{x) = fix)/ fix) 

„2 , IP -3 



S 



E-ix-^ = 0. (25) 



In the selfadjoint case, the eigenfunctions have nodes on the real axis, which leads to 
singularities in the logarithmic derivative. Thus, although the order of the equation 
can be lowered, this is not quite helpful for numerical calculations. In our case, the 
eigenfunctions corresponding to real eigenvalues are PT-symmetric (eventually up to 
an irrelevant phase factor). This implies that if s(a;; E) is the logarithmic derivative of 
an eigenfunction, then s(0; E) is purely imaginary. The large x asymptotic behaviour 
of the solution of (KS) whose real part goes to — oo as a; — > +oo is 



1 -f i 3 3 1 — i 3 

six) fa —x"^- — ^ ^Ex^^ + O (x^i ] ; x ^ +oo. (26) 

V2 4a; 2V2 l / ' 
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Now, choosing a sufficiently large X, where s(X) for the solution which decays at 
infinity is close to the asymptotic, we can integrate ( |25| ) from X backwards to x = 0. 
ff E is an eigenvalue, then 

Re[s(0; £;)] = 0. (27) 

High precision solutions of (2^) could be obtained by a simple Newton routine, starting 
at the WKB prediction (|l). 

The advantage of this approach is that we integrate a complex first order equation 
instead of the original second order one. The requirement that the imaginary part of 
s is negative allows stopping the integration if we detect a change of sign in Im(s, E). 
That energy cannot be an eigenvalue. If X is sufficiently large, the eigenvalues are 
relatively insensitive to small perturbations of the boundary condition at X. Practical 
coonsiderations require us to start with the lowest value of X, where the decaying 
solution is sufficiently close to the asymptotic behaviour (^6|) . An alternative approach 
for the large eigenvalues, would be to obtain corrections to the leading order WKB 
result 
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